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THE COMPLETE CLASSIFICATION OF FIVE-DIMENSIONAL DIRICHLET-VORONOI 
POLYHEDRA OF TRANSLATIONAL LATTICES 

MATHIEU DUTOUR SIKIRIC, ALEXEY GARBER, ACHILL SCHURMANN, AND CLARA WALDMANN 


Abstract. In this paper we report on the full classification of Dirichlet-Voronoi polyhedra and Delaunay 
subdivisions of five-dimensional translational lattices. We obtain a complete list of 110244 affine types (L- 
types) of Delaunay subdivisions and it turns out that they are all combinatorially inequivalent, giving the same 
number of combinatorial types of Dirichlet-Voronoi polyhedra. Using a refinement of corresponding secondary 
cones, we obtain 181394 contraction types. We report on details of our computer assisted enumeration, which 
we verified by three independent implementations and a topological mass formula check. 


1. Introduction 


The study of translational lattices and their Dirichlet-Voronoi polyhedra are classical subjects in crys¬ 
tallography. In 1885 Fedorov IIFed85ll (cf. IISG84I1 ) determined the five combinatorial types of possible 
Dirichlet-Voronoi polyhedra in the Euclidean 3-space These are also all the parallelohedra in that 
is, polyhedra admitting a facet-to-facet tiling of by translation. Voronoi IIVor08ll developed a theory to 
classify Dirichlet-Voronoi polyhedra for arbitrary d-dimensional Euclidean spaces His theory allows to 
classify them via a classification of Delaunay subdivisions up to affine equivalence (so called L-types). In 
this context Voronoi also came up with his famous and still unsolved conjecture, stating that every parallelo- 
hedron in is affinely equivalenf fo a Dirichlef-Voronoi polyhedron for some franslafional laffice. 

In fhis paper we reporf on fhe enumeration of fhe 5-dimensional combinatorial fypes of Dirichlef-Voronoi 
polyhedra or equivalenfly Delaunay subdivisions (Theorem l3.5l) . We find in fofal 110244 differenl combina¬ 
torial fypes and hereby go beyond fhe parfial classificafion according fo subordinafion schemes previously 
obfained by |Eng00|. In Table [3] we lisf fhe number of Delaunay subdivisions fhaf were computed so far. 
By our work, a full classification is known for d < 5 so far. Recenf parfial resulfs on primifive fypes in 
dimension 6 IIBE13II seem fo indicate fhaf a full classificafion beyond 5 dimensions is ouf of reach af fhe 
momenf. 

Our paper is organized as follows. In Secfion|2]we sfarf wifh some nofafion and background on Dirichlef- 
Voronoi and Delaunay polyfopes. Voronoi’s E-fype fheory is briefly reviewed in Section[3l We in parficular 
describe how fhe classificafion of Dirichlef-Voronoi polyhedra is reduced fo fhe classificafion of Delaunay 
subdivisions and how fhis can pracfically be done. Algorifhms and implemenfafions for our classificafion 
resulf are briefly described in Section |4] and references fo online sources are given. Additional dafa and 
fables are presented in Section [51 where we also relafe our work fo fhe fheory of confracfion fypes. 


2. Dirichlet-Voronoi and Delaunay polytopes 

Eel A denote a translational lattice in Thai is, A is a full rank-discrete subgroup of and, equiva¬ 
lenfly, can be wriffen as 

A = {Ai6i + ... + Xdbd ■ Ai,..., Aci € Z} 

wifh linearly independenf vectors bi,... ,bd £ Eafler vecfors, as well as a mafrix B wifh fhese as 
columns, are referred fo as a basis of A and we simply wrife A = HZ'^. Viewing as a Euclidean space 
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with norm | • |, the Dirichlet-Voronoi polytope of A is defined as the set of points in ! 
close to the origin than to any other element of A: 

■ad 


which are at least as 


DV(A) = |a; € : |a:| <\x — y\ for all y e a| 


2.1. General facts about polytopes. The term polytope refers to the fact that DV(A) can be described as a 
convex hull (set of all convex combinations) of finitely many points. A point that can not be omitted in such a 
description is called a vertex of the polytope. Let us briefly review some basics from the theory of polytopes 
(see IIZie951 IGriiOSII for details). A supporting hyperplane is an affine hyperplane having the property that 
the polytope is fully contained in one of the two halfspaces bounded by it. A k-dimensional face of a 
polytope is defined as a fc-dimensional intersection of the polytope with a supporting hyperplane. The 
[d — 1)-dimensional faces of a d-dimensional polytope are called/accA and vertices are the 0-dimensional 
faces. Every polytope also has a description by linear inequalities and the non-redundant ones in such a 
description are in 1-to-1-correspondence to its facets. 

Altogether, the faces of a polytope form a poset (partially ordered set, ordered by inclusion), which is 
called the/ace lattice of the polytope. Two poly topes are called combinatorially equivalent, if they possess 
the same face lattice. For instance, two 2-dimensional n-gons (which are the 2-dimensional polytopes with 
n vertices) are always combinatorially equivalent. However, the y might not be affinely equivalent, that is, 
there does not exist an affine map, mapping one to the other (see jBDSP^ldt for details on this and how to 
compute equivalence). 

We note that Engel | EngOO | uses a so called subordination scheme (sometimes called polyhedral scheme) 
which is an invariant to classify Dirichlet-Voronoi polytopes. Two combinatorially different polytopes can 
however have the same subordination scheme. In fact, several combinatorially different Dirichlet-Voronoi 
polyhedra in have the same subordination scheme. Therefore this invariant can not be used for a full 
classification of all combinatorial types. 


2.2. Affine and combinatorial types of Dirichlet-Voronoi polytopes. In dimension 2 there exist only two 
combinatorially inequivalent types of Dirichlet-Voronoi poly topes: either centrally symmetric hexagons or 
rectangles. We note that there are infinitely many affine types of Dirichlet-Voronoi polytopes. Actually, 
any centrally symmetric hexagon with vertices on a unit circle is a Dirichlet-Voronoi polytope of a lattice. 
However, they are not all affinely equivalent to each other. For instance, none of them is affinely equivalent 
to a regular hexagon (except the regular hexagon itself). We refer the interested reader for more information 
on affine types of Dirichlet-Voronoi poly topes to IIDINllirGavl4ll . 

The combinatorial types of Dirichlet-Voronoi poly topes in dimensions 3 and 4 are known as well. There 
exist five different combinatorial types of Dirichlet-Voronoi polytopes in dimension 3 and 52 different com¬ 
binatorial types in dimension 4. In this paper we report on the classification in dimension 5 and we show: 


Theorem 2.1. There are precisely 110244 combinatorially inequivalent types of Dirichlet-Voronoi polytopes 
of five-dimensional translational lattices. 


In the following we explain in more detail how to obtain the above classification result, based on Voronoi’s 
second reduction theory for positive definite quadratic forms. 

2.3. Delaunay subdivisions. The notion of Delaunay subdivisions were introduced in IIDel34ll . Here we 
give their definition and shortly describe major properties. 

Given a translational lattice A in M'^, an empty sphere S(c, r) of center c and radius r > 0 is a sphere 
such that there is no lattice point in its interior. A Delaunay cell is an intersection A D 5(c, r). A Delaunay 
polytope is a d-dimensional polytope of the form conv(A D 5(c, r)). 

The set of all Delaunay polytopes of A form a polytopal subdivision of called the Delaunay subdi¬ 
vision of A. In general, a polytopal subdivision is a non-overlapping union of polytopes that fill all of 
and such that the intersection of any two polytopes is either empty or a fc-dimensional face. DV (A) together 
with all its translates by lattice vectors form another polytopal subdivision of M'^. Both subdivisions are 
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invariant by lattice translations. The Delaunay polytopes with vertex at x G A are translates by x of some 
Delaunay polytope with vertex at 0. Thus to know the full Delaunay subdivision of a lattice A, it suffices 
to know the Delaunay polytopes with vertex 0. The centers of these Delaunay polytopes coincide with the 
vertices of DV(A). 

The Delaunay subdivision is said to be dual to the subdivision with Dirichlet-Voronoi polytopes. The 
Dirichlet-Voronoi polytope of a lattice can be obtained from the Delaunay polytopes with vertex 0 and vice 
versa: There is a bijection between the /c-dimensional faces of these Delaunay polytopes and the {d — k)- 
dimensional faces of the Dirichlet-Voronoi polytope. In particular, each d-dimensional Delaunay polytope 
corresponds to a vertex of the Dirichlet-Voronoi polytope. Moreover, the face lattice structure with respect 
to inclusion is preserved as well: If two faces of Delaunay polytopes with vertex 0 are contained in each 
other, the corresponding dual faces of the Dirichlet-Voronoi polytope are contained in each other with the 
inclusion reversed. Therefore, the classification of combinatorial types of Dirichlet-Voronoi polytopes is 
equivalent to the classification of combinatorial types of Delaunay subdivisons. 

The different combinatorial types can be derived from possible affine types. Here, two Delaunay subdi¬ 
visons, respectively lattices A and A' are affinely equivalent (are of the same affine type), if there is a matrix 
(linear map) A G GLrf(M) with A' = AK, mapping all Delaunay polytopes of A to those of A'. Note that 
two Delaunay subdivisions with different combinatorial types can not be affinely equivalent. The opposite 
could be possible though: Two different affine types of Delaunay subdivisions could possibly have the same 
combinatorial type — although we do not know of a single example among Delaunay subdivisons for trans¬ 
lational lattices at this point. In paiticular, up to dimension 5, all affine types of Delaunay subdivisons are 
not only affinely inequivalent, but also combinatorially. 


3. VORONOI’S SECOND REDUCTION THEORY 

In the following we give a short sketch of Voronoi’s second reduction theory BVorOSII . as far as it is 
necessary to describe how our classification of affine types of five-dimensional Delaunay subdivisions is 
obtained. For a more detailed description and extensions of the theory we refer to IISch09ll . 


3.1. Working with Gram matrices. The set of real symmetric positive definite matrices is denoted <S>q. 
When dealing with lattices up to orthogonal transformations, it is often convenient to work with Gram 
matrices Q = B^B G instead of using matrices of lattice bases B. Up to orthogonal transformations, 
the basis matrix B can uniquely be recovered from Q using the Cholesky decomposition. Geometrically this 
is equivalent to reconstruction of a basis knowing vector lengths and angles between them. Every positive 
definite symmetric matrix Q defines a corresponding positive definite quadratic form x Q[x\ = x^Qx on 
W^. 

In particular for studying affine fypes of Delaunay subdivisions if is convenient to use the same co¬ 
ordinates of vertices vi,... ,Vn from a fixed translational lattice A C (often A = for different 
affine images B ■ convjui,..., Vn\ of Delaunay polytopes, which we represent by a corresponding matrix 
Q G SyQ. A polytope P = conv{ui,..., Vn} with vertices Uj G A is called a Delaunay polytope of Q if 
it is d-dimensional and if there exists a center c G and a real number r such that Q[c — Vi] = for 
i = 1..., n and Q[c — v\ > for all other u G A. The set Del(A, Q) of all Delaunay polytopes of Q G 5>q 
is a polytopal subdivision of called the Delaunay subdivision of Q with respect to A. 

We speak of a Delaunay triangulation, if all the Delaunay polytopes are simplices, that is, if all of them 
have affinely independent vertices. We say that Del(A, Q) is a refinement of Del(A, Q') (and Del(A, Q') is a 
coarsening of Del(A, Q)), if every Delaunay polytope of Q is contained in a Delaunay polytope of Q'. Any 
Delaunay subdivision can be refined to a Delaunay triangulation by perturbing Q if necessary. Voronoi’s 
theory of secondary cones which we explain below gives us an explicit description of the set of positive 
definite matrices having the same Delaunay subdivision. 
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3.2. Secondary Cones and L-types. Voronoi’s second reduction theory is based on secondary cones (also 
called L-type domains) 

SC{V) = {q G 5^0 : Del(Z^g) = p} , 

which can be seen to be non-empty polyhedral cones in 5>o (which are open within their linear hull), if V 
is a Delaunay subdivision for some Q. In order to give an explicit description of SC{V) we define for an 
affinely independenf sef C C of cardinality d + 1 and a poinf w € fhe symmetric matrix 

( 1 ) Ny^w = ayVV^, 

vGV 

where the coefficients are uniquely determined by the affine dependency 

w = ^ a^v wifh 1 = ^ q;„. 
vGV uSV 

In fhe special sifuafion of C = {ui,..., being verfices of a Delaunay simplex L and w being the 

additional vertex of a Delaunay simplex L' = conv{r; 2 ,..., ru} adjacent to L, we use the notation 
Nl,L' for Nv,w In the following we use {A, B) = Trace {AB) to denote the standard inner product defined 
for fwo symmetric matrices A, B on . The following result by Voronoi gives an explicit description of a 
secondary cone in terms of linear inequalities. 

Theorem 3.1 ( llVorOSII ). Let Q be a positive definite symmetric matrix whose Delaunay subdivision D = 
Del{Z^, Q) is a triangulation. Then 

(2) SC(D) = {Q' G Q') > 0/or adj. L, L' G V}. 

This theorem of Voronoi shows that the secondary cone SC(D) of a Delaunay triangulation D is a full 
dimensional open polyhedral cone, that is, the intersection of finitely many open halfspaces. If we use weak 
inequalities > 0 in (O instead of strict inequalities, we obtain a description of the closed polyhedral cone 
SC(D). We will use these closed versions and their facial structure in the sequel. Just like for polytopes 
(cf. Section 12.11 ). faces can be defined for these closed polyhedral cones and the set of all faces forms a 
combinatorial lattice - the face-lattice of the cone. Voronoi discovered that the faces of 80(22) correspond 
to all the possible coarsenings of D. 

Two full dimensional secondary cones touch in a facet, if and only if the corresponding Delaunay trian¬ 
gulations can be transformed into each other by bistellar flips. That is we first apply a coarsening of some 
of the simplices to repartitioning polytopes (d-dimensional polytopes with d -|- 2 vertices) and then apply a 
refinement procedure. Since these changes of Delaunay triangulations are not important for what follows, 
we omit a detailed description here and refer the interested reader to IISch09ll . 

The rational closure >q of iS^q is the set of positive semidefinite quadratic forms whose kernel is 
defined by rational equations. At the core of Voronoi’s theory is the action of the matrix group GLfiZ) on 
the polyhedral tiling by closures of secondary cones: 

Theorem 3.2 (Voronoi’s Second Reduction Theory). The topological closures 80(22) give a polyhedral 
subdivision of >q on which the group GLrf(Z) acts by 80(22) C/*80(P)(7. Under this group action 

there are only finitely many inequivalent secondary cones. 

Note that one can subdivide the secondary cones into smaller cones and obtain a reduction domain for the 
action of GLfiZ) on S^q. This is the reason why Voronoi’s theory of Delaunay subdivisions and secondary 
cones is referred to as Voronoi’s second reduction theory (for positive definite quadratic forms). 

For our classification of affine types the following observation is crucial: 

Theorem 3.3. Let Q, Q' G iS^q be two positive definite matrices with Cholesky decompositions Q = B^B 
and Q' = {B'y{B') and corresponding lattices A = BZ'^ and A.' = B'Z'^. Then the Delaunay subdivisons 
of A and A' are of the same affine type if and only if Q and Q' are in GIjfiZ)-equivalent secondary cones. 

4 











Proof. We are not aware of an explicit reference for this result, so for clarity we give an argument here. First 
we note that transforming a set A and a Delaunay decomposition Del(A, Q) by a linear map A € GLd(]R) 
we get a new Delaunay decomposition Del(A', {A~^YQA~^) with vertex set A' = AK. 

Suppose now that the Delaunay decompositions of A and A' are of the same affine type. Then A ■ 
Del(A, Idd) = Del(A', {A~^YA~^) = Del(A',Idrf). Therefore 

Del(Z'^,Q) = 

= B-M-iDel(A',Idd) 

= C/Del(Z'^,Q') 

with U = B~^A~^B'. Since we have U £ GLd(Z) and therefore Q and {U~^YQ'U~^ are in 

the same secondary cone. 

On the other hand, if Q and Q' are in GLrf(Z)-equivalent secondary cones, then there exists a [7 £ 
GLrf(Z) with Del(Z'^, Q') = f7Del(Z‘^, Q). Thus 

(B')"'Del(A',Idrf) = [/B-^Del(A,Idd), 

and hence A = B'UB~^ satisfies ADel(A, Id^) = Del(A', Id^). □ 

Wifh fhe knowledge on how fo perform bisfellar flips. Theorems 13.21 and 13.31 easily lead fo an algorifhm 
fo enumerafe all affine fypes of Delaunay friangulafions in a given dimension (see Algorifhm 3 in IISch09ll ). 
For if, Schiirmann and Vallenfin developed fhe program see (secondary cone cruiser). Ifs firsf version 
from IISV05II already allowed fo reproduce fhe known classification of all GLrf(Z)-inequivalenl Delaunay 
friangulafions up fo dimension d = 5. We will use fheir resulf, respecfively fhe oufpuf of fhe program see. 

Beginning wifh dimension 6 fhe number of inequivalenf Delaunay friangulafions sfarfs fo explode. Af 
fhe momenf, we sfill do nof know how many inequivalenf friangulafions we have fo expeef in dimension 6. 
Baburin and Engel IIBE13II reporf fhaf fhey found 567.613.632 ones so far. 


3.3. Enumeration of all Delaunay subdivisions. Arbifrary Delaunay subdivisions are limiting cases of 
Delaunay friangulafions. Their secondary cones occur on fhe boundaries of full-dimensional secondary 
cones of Delaunay friangulafions. The following fheorem seems fo be folklore. One can find a proof for 
example in Proposition 2.6.1 of IIVal03ll : 

Theorem 3.4. Let B be a Delaunay triangulation. 

(1) A positive definite symmetric matrix Q lies in SC{D) if and only ifD is a refinement o/Del(Q). 

(2) If two positive definite symmetric matrices Q and Q' both lie in SC{D), then Del((5 + Q') A a 
common refinement ofBe\{Q) and Del(Q'). 

We nofe fhaf fhis fheorem can be exfended fo posifive semidefinife symmefric mafrices in fhe rational 
closure 

>0 ofSU. Eor fhose among fhem which are nof posifive definife, one can define a polyhedral 
Delaunay subdivision wifh unbounded polyhedra. Eor defails we refer fo Chapfer 4 of HSchOQII . 

By Theorem 13.41 fhe classification of all inequivalenf Delaunay subdivisions is equivalenf fo fhe classi¬ 
fication of all inequivalenf secondary cones. In order fo prove our Theorem 12.11 we show fhe following 
equivalenf resulf: 


Theorem 3.5. In dimension 5 there are 110244 affine types of Delaunay subdivisons. Equivalently, there 
are that many secondary cones of positive definite quadratic matrices in up to Glj^{'L)-equivalence. 


3.4. Related works. At this point, we should point out that there is a parallel theory that considers a single 
Delaunay polytope in a lattice, irrespective of the other Delaunay polytopes in the tessellation. This theory 
is exposed in BDEIOII and recent developments can be found in BPS 1611 . The possible Delaunay polytopes of 
dimension 5 were classified in IIKon02ll in terms of 138 combinatorial types. The classification in dimension 
6 in IIDut04ll gives 6241 combinatorial types. 

In IISch09ll (cf. Table 2 on page 60) it is reported that Engel |Eng00| found 179372 inequivalenf five¬ 
dimensional Delaunay subdivisions. This, however, is unfortunately a misinterpretation of Engel’s result 
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who classifies so called contraction types (of parallelohedra). From these contraction types, he derives 
103769 “combinatorial types”. These types are not the true combinatorial types that are classified here 
however, but a coarser notion, which classifies parallelohedra in dimension 5, or equivalently Delaunay 
subdivisions, up to their subordination schemes. The subordination scheme of a d-dimensional polytope P 
is a list of numbers, containing for every k = 2,... ,d — 1 and for every n, the number of {k — l)-faces 
of P incident to exactly n of the /c-faces of P (see Section 4 of |Eng00 1 for details). Thus, the subordination 
scheme encodes certain properties of the face lattice of a polytope, but not the whole face lattice. Two 
combinatorially different polytopes can have the same subordination scheme. They may even be the same for 
different affine types of Dirichlet-Voronoi polytopes, having even secondary cones of different dimension. 
In fact, during our work we discovered two such examples for d = 5. 

Note that combinatorial types of polytopes can only truely be distinguished, by checking whether or not 
their face lattices are different. It has been shown in IIKS03II that the incidence relations between vertices 
and facets of two polytopes are sufficient to distinguish their face lattices. Practically such differences can 
be checked using graph isomorphism software as we describe in the next section. Invariants like the number 
of faces of a given dimension or the subordination scheme used by Engel may be useful in computations, 
for instance when limiting the number of equivalence tests. However, such invariants are not sufficient 
for complete enumerations. Engel’s invariant appears to distinguish the known 52 combinatorial types in 
dimension 4, but it does not distinguish types in any dimension greater or equal to 5. While it is conceivable 
that the subordination scheme could be extended to better distinguish between types, it should never be used 
alone without checking for equivalence since there is always the possibility that non-isomoiphic structures 
have the same invariant. 


4. Algorithms and Implementations 

Before we explain the details of our computations for d = 5, we start with some general observations, 
which are valid in all dimensions and quite useful for practical purposes. 

4.1. Using reduced generators and central forms. Each secondary cone, respectively its closure is given 
by a finite list of linear inequalities (coming from Voronoi’s regulators, cf. Theorem 13.11) . Erom it one 
can obtain a number of generating rays. In fact, one of these descriptions (by rays or inequalities) can be 
obtained from the other by a polyhedral representation conversion. Since all of the involved inequalities 
involve rational numbers only, we may assume that the generators for rays are given by integral vectors 
(matrices in <S'^), with coordinates having a greatest common divisor (gcd) of 1. We refer to these generators 
as reduced (or normalized) generators. As we are using Theorem 13.41 for the classification of Delaunay 
subdivisons, we only need to consider closures of secondary cones which are faces of closures of full¬ 
dimensional secondary cones. All such faces are themselves generated by a subset of the reduced generators 
of the full dimensional cone. 

Having reduced generators iii,..., of a secondary cone SC (respectively its closure), we define a 
central reduced (or normalized) form of the secondary cone as the sum Q{SC) = H to 

that two secondary cones SC and SC' are GLrf(Z)-equivalent if and only if Q(SC) and Q{SC') are GLrf(Z)- 
equivalent. Hence, for the classification of secondary cones up to GLrf(Z)-equivalence we can equally well 
classify their central reduced forms up to GLrf(Z)-equivalence. 

4.2. Testing equivalence of forms and use of invariants. Testing GLrf(Z)-equivalence of central reduced 
forms can be done with the Plesken-Souvignier algorithm IIPS97II . Their initial implementation is available 
at IPS 95 II and is pait of computer algebra software such as IIMAG061IGAP151 . The algorithm works by 
building a finite set of vectors that is canonically defined by a given positive definite matrix and spans Z'^ as 
a lattice. Eor a given norm bound n and a positive definite matrix Q let 

S{Q, n) = jn G s.t. Q[v] < n| . 
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Then we take the smallest n such that S{Q,n) spans as a lattice and call the vector set Can((5). 

As testing GLd(Z)-equivalence of central reduced forms is computationally quite involved, one needs to 
reduce the number of such tests as much as possible since the final number of forms is M = 110244 and so 
the total number of isomorphism tests is a priori M{M — l)/2. The basic idea is to use invariants to reduce 
the number of tests. Some invariants come naturally from the form (5(SC) such as its determinant and size 
of Can((5(SC)). Other possible invariants are related to the secondary cone SC under consideration. For 
example the dimension of SC or its number of generating forms Ri, ..., R^. Further invariants are the rank 
of Rk and so on. Rather surprisingly, the most efficient invariant tends to be the determinant of Q(SC). 

4.3. Putting it all together for five dimensions. Now, finally, lef us pul Ihe pieces above logelher, lo 
describe Ihe algorilhm behind our classification resull for d = 5. To show Theorem 13.51 wilh computer 
assislance, we can use Voronoi’s Iheory. We slait from Ihe secondary cones of Ihe 222 known Delaunay 
Iriangulalions. Those were classified in IIBR731IRB781 bul Ihe classification was incorrecl and a final correcf 
classification was oblained in IIEG02II which we have independenlly confirmed in IISV061IDSG091 . These 
open polyhedral cones are full dimensional in and Iherefore have dimension 15. Their closure is given 
by a lisl of non-redundanl linear inequalities. From Ihis lisl, we can oblain Ihe reduced generafors of each 
cone and also a description by generafors and by equalions / inequalilies for each of Iheir facels. These 
facefs are Ihemselves closures of 14-dimensional secondary cones which correspond fo Delaunay subdi¬ 
visions lhal are a Irue coarsening of Ihe considered Delaunay Iriangulafion al hand. Some of Ihem may be 
GLrf(Z)-equivalenl, so for our classification, we have fo oblain a lisl of GLrf(Z)-inequivalenl 14-dimensional 
secondary cones in <S>q from Ihem, using Iheir cenlral reduced forms. In a nexl step, we oblain a lisl of 
GLrf(Z)-inequivalenl 13-dimensional secondary cones from our lisl of 14-dimensional secondary cones in a 
similar way. We continue Ibis process until we subsequenlly oblain a full lisl of GLrf(Z)-inequivalenl cones 
of dimensions 15,..., 1. See Tabled] for the number of secondary cones oblained in each dimension in Ibis 
way. 


4.4. Practical Implementations. The computer code of our first implementation in Haskell of the al¬ 
gorithm described above, together with a detailed documentation (in German) is available at the web¬ 
page HWallSbl . In particular, data of the full classification can be obtained at HWallSal . with a matrix 
of a central reduced form for each secondary cone in 5 ^q. 

Our second implementation used the GAP package polyhedral UPS 1511 with some external calls to 
isom BPS 95 II for equivalence tests and Irs BAvil5ll for polyhedral representation conversions. In our third 
implementation, we adapted the program see. In its latest version IIGSV15II we included the program isom 
to produce all secondary cones of a given dimension. 

In order to avoid the dependency on isom in all three implementations, we also performed equivalence 
computations with nauty BMcKlTl . applied to test equivalence of the sets Can(Q(SC)) of vectors, by 
using the method explained in Section 3.4 of llBDSP'*~14l . Overall, the full computation, its resulting data 
and in particular the numbers in Table [T] were all sufficiently well cross-checked. All calculations yield 
the same results and due to the different nature of our three programs we can be certain of the obtained 
classification, although the computations are large and quite involved. 

We can use the obtained results for a computational proof of our main Theorem 12.11 by showing that 
all Delaunay subdivisions, respectively the corresponding Dirichlet-Voronoi polytopes, are combinatorially 
inequivalent. This is shown by checking if their face-lattices are non-isomorphic. Since the face-lattice of a 
polytope is determined by the incidence graph of vertices and facets, we can check if these graphs are non¬ 
isomorphic. These isomorphism checks can be performed using for instance graph isomorphism software 
like nauty BMcKlTl . We computed “canonical forms” for each of the graphs with nauty and then used 
mdSsum (a special hash function) for each of them in order to decide computationally (in a reasonable 
amount of time) that they are all different. 
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5. Tables and Data 


We provide the following tables, containing additional information: Table [Ugives the number of inequiv¬ 
alent secondary cones by their dimension. Table [2]gives the number of secondary cones by their number of 
rank 1, 4 or 5 extreme rays. Table|3]gives the known numbers of inequivalent secondary cones (all combina¬ 
torial types) and full-dimensional secondary cones (primitive types), together with a reference where these 
results can be found. Table |4] gives the number of secondary cones according to their dimension and their 
number of extreme rays. Table Ogives the number of secondary cones that cannot be extended to a higher 
dimensional cone by a pyramid construction with a rank-1 extreme ray. Table [ 6 ] gives the frequencies of 
occuning Bravais groups according to the nomenclature of UCAROSl . Table |7] and [ 8 ]relate our classification 
to notions in the theory of contraction types as developed in | EngOO[ . In the following we provide some 
background information (see also IIDSGM14II ). 


5.1. Fundamental faces and irreducible cones. For a given secondary cone SC with generating rays 
..., we define the fundamental face F(SC) fo be fhe smallesf face of SC fhaf confains all fhe gen- 
erafors Ri of rank greater fhan 1 . The face F{SC) may be reduced fo zero in which case SC is generated 
by rank-1 mafrices only. From IIFR94II we know fhaf fhe number of generators is equal fo fhe dimension 
of fhe secondary cone in fhis case and fhaf fhis case is equivalenf fo fhe Dirichlef-Voronoi polyfope being 
a zonofope and fo fhe Delaunay subdivision being fhe connecfed region of a hyperplane arrangemenf. Up 
fo GL 5 (Z)-equivalence, we found 81 secondary cones of fhis kind, corresponding fo differenf zonofopes in 
dimension 5. 

If F(SC) is nonfrivial (non-zero) fhen fhe sfrucfure of fhe secondary cone is more complex. For a sec¬ 
ondary cone SC we have a decomposifion of fhe form 

h 

SC = F(SC) + ^M+p(r;i), 

i=l 

wifh p{vi) = Viv\ fhe rank-1 matrix (form) associated to a vector Vi. Our computations show that we have 
dim SC = dim F(SC) -I- h which means that SC is obtained by a sequence of h pyramid constructions over 
F(SC). By a pyramid construction we mean an extension to a higher dimensional secondary cone by adding 
a rank -1 generating ray. 

If F(SC) does not contain any positive definite matrices (and hence lies in the boundary of iS|,q), then in 
dimension 5 there is only one possibility: F(SC) has only one extreme ray that corresponds to the D 4 root 
lattice, which we denote by FD 4 . Up to GL 5 (Z)-equivalence, we found 424 different combinatorial types 
of secondary cones of the form '^+p{vi)- Note that ^04 itself is not a secondary cone, since it 

does not contain any positive definite forms. By our computation, all such cones have their dimension equal 
to their number of generators. 

The fundamental cones F(SC) may themselves contain rank 1-forms. For example, there exist two 
secondary cones of dimension 3 with 4 generators each, 3 of rank 4 and one of rank 1 (see Section 5 of 
IIDSHS15II ). If F(SC) contains only forms of rank higher than 1 then according to the terminology of 
HFngOO I it is totally zone contracted. If a secondary cone satisfies SC = F(SC) fhen if is called irreducible. 
Table |7] and [ 8 ] give key information on irreducible secondary cones we found. 


5.2. Contraction types. In [EngOOl the notion of a contraction type is introduced. This notion is distinct 
from secondary cones and gives a further refinement of them. That is, if we have a secondary cone SC that 
is irreducible but not totally zone-contracted and has rank -1 forms pi, ..., pm, then we can decompose it 
into a number of contraction cones (also called contraction domains) SCj -|- ^+Pj with SCj a totally 
zone-contracted secondary cone. For example the 3-dimensional cone SC with symbol L^L^pi in Table [ 8 ]is 
a cone over a square (combinatorially) with vertices corresponding to pi, Li, L 3 and Li. We can decompose 
it into two isomorphic 3-dimensional cones (over triangles) of the form L 1 L 3 -I-M+P 1 and one 2-dimensional 
cone of the form L 3 -|- M+pi. 

















For other cones the decomposition can be more complicated. Given an irreducible secondary cone SC, 
let Ri be the cone of its extreme rays of rank 1. We define S to be the set of all totally zone contracted 
irreducible cones whose rays are also rays of SC (of rank greater than 1). Then our computation shows that 
SC can be decomposed into contraction cones S + Ri with 5 G 5. 

The decomposition of an irreducible secondary cone SC into contraction cones, induces a decomposition 
of any secondary cone obtained by adding rank-1 forms. Overall, we thus obtain a decomposition into 
contraction cones that is finer than the decomposition by secondary cones. For secondary cones SC whose 
fundamental face F(SC) is totally zone-contracted there is no difference. But for other irreducible secondary 
cones the contraction types form a strictly finer decomposition. The total number of contraction types 
that we obtain is 181394. The number of contraction cones by their dimension is given in Table [U In 
Table [8] we give for each irreducible secondary cone D the number of types of contraction cones contained 
in L) -|- M.^p{vk)- We note that in [EngOOl the number of contraction cones is reported to be 179372. 

This discrepancy is most likely due to the different notion of equivalence via “subordination schemes” used 
there. 


5.3. Euler Poincare characteristic check. Another key check of the correctness of our enumeration is to 
use the Euler Poincare characteristic. We have the formula 




Stab(F)| 


= 0 


where the sum is over the representatives of cones with respect to the action of GL„(Z). This kind of 
formula co mes from the Euler Poincare characteristic of discrete groups, i.e. x(GL„(Z)) = 0 for n > 3. 
See IIBro94[|DSGG~*~l^ for more details. 

Both, our enumeration of secondary cones and our enumeration of contraction cones satisfy this con¬ 
dition, which is yet another strong indication of the correctness of our enumeration. Eor example for the 
secondary cones, if we regroup the cones by their dimension, this gives us the following non-trivial identity: 

293 I 7463 _ 9^ i 56927 _ 5146751 

I 8329297 ^’^^ °3341 I l&i)783 ^’^ ?03083A^ ^f° 13879537 

fIF&53 I lll'^727 _ ^1^95 , 48^1)^ 8923 

180 d- 288 288 d- gg 144 u. 

This kind of mass formula provides a highly non-trivial check of the correctness of an enumeration as any 
error on a single entry or on a single stabilizer would turn the formula wrong. 
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Table 1. Number of GL 5 (Z)-inequivalent secondary cones and contraction cones in <S|.q 
by their dimension. 


n 

nr. sec. c. 

nr. cont. c. 

n 

nr. sec. c. 

nr. cont. c. 

1 

7 

7 

9 

21132 

33085 

2 

37 

39 

10 

22221 

37601 

3 

146 

161 

11 

18033 

32821 

4 

535 

613 

12 

10886 

21292 

5 

1681 

2021 

13 

4713 

9709 

6 

4366 

5543 

14 

1318 

2787 

7 

8 

9255 

15692 

12512 

22806 

15 

222 

397 


Table 2. Number of GL 5 (Z)-inequivalent secondary cones in <S>g by number of rank-A: 
generating rays. In line i, the rank-k-column, k = 1,4,5, contains the number of secondary 
cones which have i generating rays of rank k. (There exist no generating rays for k = 2, 3.) 

# Generating rays 


(of particular rank) 

rank-1 

rank-4 

rank-5 

0 

82 

51900 

1572 

1 

410 

35316 

15421 

2 

1658 

21574 

32939 

3 

5029 

1354 

26811 

4 

11301 

0 

19302 

5 

18923 

100 

6841 

6 

23802 

0 

3662 

7 

22411 

0 

2150 

8 

15528 

0 

950 

9 

7744 

0 

285 

10 

2699 

0 

170 

11 

548 

0 

38 

12 

97 

0 

76 

13 

9 

0 

0 

14 

2 

0 

0 

15 

1 

0 

9 

16 

0 

0 

18 


12 












Table 3. Number of primitive and all combinatorial types of Delaunay subdivisions, re¬ 
spectively corresponding GL„(Z) -inequivalent secondary cones. 


n 

Primitive types 

All combinatorial types 

2 

1 

2 

3 

1 IIFed85ll 

5 IIFed85ll 

4 

.3 IIVor08ll 

52 IIDel29llSto751 

5 

222 IIBR73II 
IIRR781IRG02I 

110244 

6 

> 567.613.632 mF4 3ll 



Table 4. Number of secondary cones according to dimension (at most 15) and number of 
generators (at most 26). 


1 

2 

3 

4 

5 

6 
7 


37 


9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


3 4 


144 

2 517 

17 


1 


5 6 7 8 9 10 11 12 13 14 15 


1595 


81 

4041 




1 

301 

8266 




12 

887 

13354 


3 


62 

2007 

16862 

1 

11 

2 

222 

3461 


1 

36 

13 

557 



2 

89 

50 




7 

182 


19 

1 


16358 


4443 

11989 





944 

4259 

6395 




122 

1103 

2945 

2346 



305 

181 

857 

1449 

526 


43 

403 

173 

430 

456 

62 


80 

390 

102 

120 

84 

5 


92 

274 

35 

13 


15 


72 

122 

5 



30 


29 

33 




34 


13 

1 

3 



23 

6 


4 

6 

7 


6 


13 























Table 5. Number of GL5(Z)-inequivalent secondary cones in which are not extend¬ 
able to a higher dimensional secondary cone by adding a rank-1 generating ray. 

Dimension_10 11 12 13 14 15 

# Secondary cones 1 12 40 142 266 222 

Table 6 . Frequency of occurence of Bravais groups, “name” is the standard name from 
the GAP package IICAR08L “order” is the size of the point group of corresponding lattices, 
“frequency” is the number of secondary cones that are symmetric with respect to the group. 


name 

order 

frequency 

name 

order 

frequency 

name 

order 

frequency 

1 ,1,1,1,1 :1 

2 

105301 

i,i;i;i;i :17 

16 

5 

4-l;l :2 

768 

1 

1 ,1,1,1;1 :2 

4 

4155 

3;1,1 :2 

96 

4 

4-l;l :3 

2304 

1 

1 ,1,1;1;1 :6 

8 

159 

3;1,1 :5 

96 

4 

5-1 :3 

3840 

1 

2 -2;l,l,l :2 

12 

137 

2 -l;l,l;l :6 

32 

4 

5-2 :3 

1440 

1 

1 ,1,1;1,1 :2 

4 

112 

l;i;i;i;i :8 

32 

4 

3;1;1 :4 

192 

1 

1,1,1;1;1 :4 

8 

90 

1 ,1,1;1,1 :1 

4 

3 

4-1;1 :4 

768 

1 

1,1,1;1;1 :5 

8 

39 

i,i,i;i;i :i 

8 

3 

2-2;2-2;l :5 

72 

1 

1 ,1,1,1;1 :1 

4 

34 

2-2;2-2;l :3 

72 

3 

2 -l;l;l;l :6 

64 

1 

2 -l;l,l,l :2 

16 

31 

1 ,1;1;1;1 :10 

16 

3 

2-l;l;l;l :7 

64 

1 

2 -2;l,l;l :6 

24 

31 

4-3;l :3 

240 

2 

2-2;l;l;l :7 

48 

1 

1,1;1;1;1 :15 

16 

20 

2-2;l,l;l :4 

24 

2 

3;1;1 :7 

192 

1 

i,i;i,i;i :3 

8 

14 

i;i;i;i;i :5 

32 

2 

2 -l;l;l;l :8 

64 

1 

1,1;1;1;1 :13 

16 

12 

2-2;l,l;l :5 

24 

2 

2 -l;l;l;l :11 

64 

1 

3;1,1 :3 

48 

10 

3;1;1 :12 

192 

2 

i;i;i;i;i :12 

32 

1 

1 ,1;1;1;1 :6 

16 

8 

1;1;1;1;1 :13 

32 

2 

2 -l;l;l;l :12 

64 

1 

3;1;1 :8 

96 

7 

1 ,1;1,1;1 :1 

8 

1 

1;1;1;1;1 :15 

32 

1 

1 ,1,1;1;1 :2 

8 

6 

1 ,1;1;1;1 :i 

16 

1 

1 ;1;1;1;1 :16 

32 

1 

2-l;l,l;l :4 

32 

6 

l;i;i;i;i :i 

32 

1 




1 ,1;1,1;1 :6 

8 

6 

3;1;1 :2 

192 

1 
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Table 7. Information about the 82 totally zone-contracted secondary cones, “dim” is the 
dimension of the secondary cone SC, “generator” gives the type of the extreme rays, “sym¬ 
bol” gives the number of facets and vertices of the corresponding Dirichlet-Voronoi poly¬ 
topes and “nb sec. c.” gives the number of secondary cones having SC as their fundamental 
face. 


dim 

generator 

symbol 

nb sec. c. 

dim 

generator 

symbol 

nb sec. c. 

dim 

generator 

symbol 

nb sec. c. 

1 

Li 

40,42 

450 

3 

LiDi 

48,242 

2738 

4 


42,204 

665 

1 

L2 

42,96 

111 

3 

L2OI 

42,168 

2047 

4 

L1L2DI 

48,282 

3988 

1 

^3 

48,180 

670 

3 

isDl 

52,344 

1344 

4 

L1L3DI 

52,352 

2272 

1 

u 

50,192 

112 

3 

L7DI 

56,462 

484 

4 

L2L3DI 

52,384 

1074 

1 

L 5 

50,282 

352 

3 

L1L2D4 

48,242 

5029 

4 

L 3 LrDl 

56,470 

1160 

1 

Lq 

54,342 

324 

3 

L1L3D4 

48,254 

2436 

4 

ilL2L3D4 

52,354 

4100 

1 

Lr 

54,366 

220 

3 

L1L5D4 

50,328 

650 

4 

^1^2^504 

54,418 

1256 

2 


42,132 

1067 

3 

L2L3D4 

52,346 

2344 

4 


54,418 

1088 

2 

L1D4 

40,122 

1814 

3 

L2L5D4 

54,402 

650 

4 


50,342 

696 

2 

-L2D4 

42,132 

1825 

3 

L3L5D4 

50,334 

553 

4 

L2L3L5D4 

54,424 

1092 

2 

L3D4 

48,246 

1428 

3 

L3L7D4 

54,410 

1160 

4 

L1L2L3L3 

54,406 

1392 

2 

L5D4 

50,312 

352 

3 

L1L2L3 

52,316 

2773 

4 

L4L2L3LQ 

54,428 

856 

2 

L7O4, 

54,402 

484 

3 

L1L2L3 

54,392 

1256 

4 

L\L2L3L& 

54,438 

928 

2 

L1L2 

48,202 

2385 

3 

L1L2L6 

54,400 

758 

4 

LIL3L4L5 

50,360 

696 

2 

L1L3 

48,188 

1058 

3 

LiLjLr 

54,382 

456 

4 

L1L3L4LQ 

54,416 

786 

2 

L1L4 

50,232 

333 

3 

L1L3L4 

50,288 

516 

4 

LIL3L3LQ 

54,418 

800 

2 

L1L3 

50,298 

650 

3 

L1L3L3 

50,312 

696 

4 

LIL4L3LQ 

54,426 

928 

2 

LiLe 

54,366 

758 

3 

L1L3L6 

54,394 

856 

4 

L2L3L3L& 

54,444 

628 

2 

L2L3 

52,308 

1638 

3 

L1L4L3 

50,346 

630 

4 

L3l/4L'5t^6 

54,432 

628 

2 

L2L3 

54,376 

650 

3 

L1L4LQ 

54,388 

734 

5 

LlDl 

42,240 

100 

2 

L2L6 

54,376 

324 

3 

LiL^Lq 

54,404 

928 

5 

LiLlDl 

48,322 

689 

2 

L3L4 

50,280 

318 

3 

L2L3L3 

54,398 

1092 

5 

L1L2L3DI 

52,392 

1815 

2 

L3L3 

50,304 

553 

3 

L2L3LQ 

54,420 

582 

5 

LiL^LrDl 

56,478 

1088 

2 

LsLq 

54,386 

582 

3 

L2L5LQ 

54,422 

553 

5 

L4L2L3L3D4 

54,432 

1392 

2 

L3L7 

54,374 

490 

3 

L3L4L3 

50,352 

553 

5 

L1L2L3L3LQ 

54,452 

800 

2 

L4L5 

50,330 

348 

3 

L3L4L6 

54,408 

531 

5 

L1L3L4L3LQ 

54,440 

800 

2 

L4L6 

54,364 

318 

3 


54,410 

628 





2 

LsLe 

54,388 

553 

3 

L4L5LQ 

54,410 

553 
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Table 8. Information about the 125 inequivalent irreducible secondary cones, which are 
not totally zone-contracted. Same convention as in Table |7j in addition pi denotes an ex¬ 
treme ray of rank 1 and “nb cont. d.” is the number of contraction cones corresponding to 
this irreducible component. 


dim 

generator 

symbol 

nb sec. c. 

nb cont. c. 

dim 

generator 

symbol 

nb sec. c. 

nb cont. c. 

3 

LfLspi 

48,196 

566 

2047 

7 

L 3 L 4 Lip\ 

58,536 

27 

72 

4 

LiL^L^pI 

50,320 

205 

3988 

7 

L\LIL 3 D 4 pI 

50,378 

73 

3030 

4 

LfL 3 D 4 pi 

48,262 

1240 

1074 

7 

LiL 2 LlL 3 pi 

54,442 

134 

639 

4 


54,390 

174 

665 

7 

LILIL4L5PI 

50,396 

73 

1274 

4 

L1L2L3P1 

52,324 

1423 

1092 

7 

LiqLsLepl 

54,454 

33 

820 

4 

LIL3L4P1 

50,296 

274 

1256 

7 

L\L 2 L 3 L 3 Lep\ 

54,474 

164 

605 

4 

LlLsLspi 

50,320 

205 

615 

7 

lIlIlIupI 

54,464 

74 

1000 

4 

L4L3L6P1 

54,402 

358 

4100 

7 

LiLiqi^pi 

54,464 

148 

740 

4 

LiL^Lspi 

50,326 

182 

3503 

7 

LlL 3 L 4 L 3 LePi 

54,462 

150 

207 

4 

LsL^Lspi 

54,434 

203 

3999 

7 

LiI/2L3LgL6Pi 

54,484 

121 

814 

5 

LiLsLqpI 

54,412 

97 

615 

7 

LiLiLfLep? 

54,480 

22 

261 

5 

L1L3L5D4PI 

50,350 

205 

1188 

7 

LiL 3 L 4 L‘^LqP^ 

54,472 

121 

1036 

5 

LlLlL^pl 

50,334 

298 

5895 

7 

LIL2LIL3D4PI 

54,454 

606 

153 

5 

L1L2L3LSP1 

54,414 

396 

492 

7 

L\L 2 L 1 L 3 Lqp\ 

54,474 

200 

750 

5 

L1L3L4L5P1 

50,368 

197 

492 

7 

LIL 2 L 3 LIL< 3 pI 

54,484 

34 

605 

5 

LiLiLzLepI 

54,426 

164 

689 

7 

LfL|L4L5L6P? 

54,462 

200 

1000 

5 

LiLsL^Lepi 

54,432 

164 

1815 

7 

LIL3L4LIL6PI 

54,472 

34 

740 

5 

LILsDIpi 

52,360 

1168 

3279 

7 

LiL 2 LlLlLepl 

54,490 

148 

207 

5 

LfL|L7D4pi 

54,426 

396 

100 

7 

LiL^L/4^L‘^Lqp^ 

54,478 

148 

639 

5 

L1L2L3D4P1 

52,362 

2060 

1392 

8 

LlLlLlL^pl 

54,478 

34 

320 

5 

L^LsLsDipi 

50,350 

205 

553 

8 

lIlIlIupI 

54,478 

47 

1274 

5 

L1L3L5D4P1 

50,356 

182 

1092 

8 

LlLlLlLepf 

54,488 

38 

337 

5 

L 1 L 2 I/ 3 L 5 P 1 

54,414 

396 

958 

8 

LlLlLpept 

54,478 

43 

814 

5 

L\L2L2,LqPI 

54,436 

358 

480 

8 

L1L2LIL5L6PJ 

54,488 

64 

487 

5 

LIL3L4LSP1 

50,368 

205 

1490 

8 

LiL'^L 4 L 5 Lepf 

54,476 

57 

285 

5 

L 7 L 3 I/ 4 L 6 P 1 

54,424 

327 

990 

8 

LiL 3 L 4 Llpf 

58,544 

28 

77 

5 

LIL3L5L6P1 

54,426 

228 

291 

8 

LIL2LIL5D4PI 

54,468 

134 

261 

5 

LiL 2 L'^L 5 Pi 

54,420 

352 

546 

8 


54,488 

33 

1036 

5 

LiL'^LiLs.pi 

50,374 

182 

800 

8 

L\L^L 4 L 3 Lep\ 

54,476 

33 

753 

5 

LiL'^L^Lepi 

54,432 

128 

628 

8 

LfL2L§L|L6Pi 

54,498 

74 

153 

5 

LiLsLlLepi 

54,442 

178 

328 

8 

LIL2LIlIL6pI 

54,498 

148 

575 

5 

L 2 L 3 L'^Lepi 

54,468 

203 

474 

8 

LfL|L4L|L6Pi 

54,486 

74 

814 

5 

L 3 L 4 LIL 6 P 1 

54,456 

203 

591 

8 

LIlIL 4 LIL 6 pI 

54,486 

148 

261 

6 

LlLepf 

54,430 

34 

92 

8 

L 1 L 2 L 3 LIL 6 P 1 

54,514 

22 

1036 

6 

LMLspI 

50,348 

73 

1188 

8 

ij 1 X/g 7^4 Z/g 

54,502 

22 

153 

6 

LILzL^LqpI 

54,440 

164 

492 

9 

LiL 4 L 3 L 3 p\ 

54,502 

16 

337 

6 

LiL 2 L 5 Lepl 

54,446 

97 

492 

9 

Lp^LlLepl 

54,502 

38 

487 

6 

LiLsLILqpI 

54,450 

121 

2619 

9 

LlL 2 LlLlLepl 

54,512 

34 

48 

6 

L1L4L5L6P1 

54,434 

93 

1092 

9 

LlLlL 4 LlLfipl 

54,500 

30 

753 

6 

LlLlL^Dipl 

50,364 

298 

958 

9 

LIL3L4LIPI 

58,552 

11 

575 

6 

L1L2L3L5D4PI 

54,440 

396 

1490 

9 

LlL2LlLlLep\ 

54,512 

47 

905 

6 

LIL2LIL3PI 

54,428 

606 

3030 

9 

L\LlL 4 LlLepi 

54,500 

47 

300 

6 

LiqiiLspl 

50,382 

298 

639 

9 

LlL2LlLlLep\ 

54,522 

38 

68 

6 

LIlILsLqpI 

54,440 

200 

291 

9 

LlL2LlLlLepf 

54,512 

43 

487 

6 

LlLsLlLepl 

54,450 

34 

820 

9 

LlLlL 4 LlL&p\ 

54,510 

38 

753 

6 

L1L2L3L5L6PI 

54,460 

164 

605 

9 

LlLlL 4 LlLep\ 

54,500 

43 

575 

6 

LiL 2 L 3 L 5 Lepl 

54,466 

164 

628 

10 

L 3 L 4 L 6 P 1 

54,452 

6 

18 

6 

LiL^L^Lqp^ 

54,456 

148 

328 

10 

LtL|L|L6pf 

54,526 

9 

70 

6 

LIL3L4L3L6PI 

54,448 

164 

1000 

10 

LiL 3 L 4 L 3 L‘^p\ 

58,582 

14 

905 

6 

L 1 I/ 3 L 4 L 5 Z/ 6 P 1 

54,454 

150 

474 

10 

LfL2L|LiL6pf 

54,536 

38 

186 

6 

LfLlLrDlpi 

56,486 

396 

740 

10 

LfL|L4Lii6pf 

54,524 

38 

905 

6 

L\L 2 L 3 D 1 pi 

52,400 

933 

207 

11 

LaLIpI 

50,468 

3 

30 

6 

L1L2L3L3D4P1 

54,440 

396 

492 

11 

LiL3L4L3Lgpf 

54,524 

8 

40 

6 

L 1 L 2 L 3 L 5 D 4 P 1 

54,446 

352 

450 

11 

L|L4Lgp| 

58,580 

6 

no 

6 

L\L 2 L 3 L 3 Lepi 

54,460 

228 

2420 

11 

LIL3L4L3LIPI 

58,590 

10 

7 

6 

L^L 3 L 4 L 3 Lepi 

54,448 

228 

279 

11 

L\L2LlLlL(ip\ 

54,560 

9 

186 

6 

LiL 2 L‘^L 3 Lepi 

54,466 

128 

1490 

11 

L\LlL 4 LlL&p\ 

54,548 

9 

186 

6 

LiL 2 L 3 L^Lepi 

54,476 

178 

628 1 

1 12 

L1L4LIL3PI 

54,548 

4 

49 

6 

LiL'^L 4 L 3 Lepi 

54,454 

128 

328 

12 

7j 1 7^3 7j 4X/5 TjgP j 

58,604 

7 

20 

6 

LiL 3 L 4 L^Lqpi 

54,464 

178 

474 

13 

LlLlL 4 LlLlp\ 

58,628 

4 

55 

7 

LiLlLsLepf 

54,454 

64 

92 

13 

LIL3L4LILIPI 

58,628 

3 

27 

7 

L 2 LlL 3 p\ 

54,464 

34 

320 

15 

L|L4L|p}2 

62,708 

1 

4 

7 

L^L^lLqp-^ 

54,452 

27 

72 



















